Effect of quantum fluctuations on even-odd energy difference in a Cooper-pair box. 
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We study the effect of quantum charge fluctuations on the discrete spectrum of charge states of a 
small superconducting island (Cooper-pair box) connected to a large finite-size superconductor by 
a tunnel junction. In particular, we calculate the reduction of the even-odd energy difference 5E 
due to virtual tunneling of electrons across the junction. We show that the renormalization effects 
are important for understanding the quasiparticle "poisoning" effect because SE determines the 
activation energy of a trapped quasiparticle in the Cooper-pair box. We find that renormalization 
of the activation energy depends on the dimensionless normal-state conductance of the junction g T , 
and becomes strong at g T 3> 1. 



o 
o _ 

u • 

Of 

3 ; 
+J . 



i 

-a 
c 

o 
o 



(N 
> 

(N 

o 

-I— > 



i 

£ 

O 

o 



X 



Recently, superconducting quantum circuits have at- 
tracted considerable interest (see [l||2| and references 
therein). From the viewpoint of quantum many-body 
phenomena, these circuits are good systems to study the 
effect of quantum fluctuations of an environment on the 
discrete spectrum of charge state a 3 i 4|5 ' 6 i 7 ' (similar to the 
Lamb shift in a hydrogen atom) . While most of the stud- 
ies of superconducting nanostructures focus on smearing 
of the charge steps in the Coulomb staircase measure- 
ments 8 , here we consider another observable quantity - 
even-odd-electron energy difference SE in the Cooper- 
pair box (CPB). This quantity is important for under- 
standing the quasiparticle "poisoning" effect ^ 10 ' 11 ! 12 ' 13 , 
and it has been recently studied experimentall y 14 ! 15 . It 
was conjectured that SE may be reduced in the strong 
tunneling regime g T = R q / Rn > 1 by quantum fluctua- 
tions of the charge^ 4 -. Here R q and Rn are the resistance 
quantum, R q = h/e 2 , and normal-state resistance of the 
tunnel junction, respectively. 

In this paper, we study the renormalization of the dis- 
crete spectrum of charge states of the Cooper-pair box by 
quantum charge fluctuations. We show that virtual tun- 
neling of electrons across the tunnel junction may lead to 
a substantial reduction of the even-odd energy difference 
SE. We consider here the case of the tunnel junction 
with a large number of low transparency channels 1 ^. 

The dynamics of the system is described by the Hamil- 
tonian 



H = R„ 



Hbcs 



BCS 



(1) 



Here H^ cs and -ff^cs are BCS Hamiltonians for the CPB 
and superconducting reservoir; H c — E c (Q/e — N g ) 2 
with E c , N g and Q being the charging energy, dimen- 
sionless gate voltage and charge of the CPB, respectively. 
The tunneling Hamiltonian H T is defined in the conven- 
tional way. We assume that the island and reservoir are 
isolated from the rest of the circuit; i.e. total number 
of electrons in the system is fixed. At low temperature 
T < T* , thermal quasiparticles are frozen out. (Here 
T* = ln (A/<5) w ith A and 8 being superconducting gap 
and mean level spacing in the reservoir, respectively). 
If total number of electrons in the system is even, then 
the only relevant degree of freedom at low energies is 



the phase difference across the junction ip. In the case 
of an odd number of electrons a quasiparticle resides in 
the system even at zero temperature. The presence of 
le-charged carriers changes the periodicity of the CPB 
energy spectrum (see Fig. [T]) since an unpaired electron 
can reside in the island or in the reservoir. Note that 
at N g — 1, a working point for the charge qubit, the 
odd-electron state of the CPB may be more favorable re- 
sulting in trapping of a quasiparticle in the islan d 14 ^ 15 ' 17 . 
In order to understand energetics of this trapping phe- 
nomenon, one has to look at the ground state energy dif- 
ference SE between the even-charge state (no quasipar- 
ticles in the CPB), and odd-charge state (with a quasi- 
particle in the CPB): 



SE = E cvcn {Ng = l)-E odd {N g = l), 



(2) 



see also Fig. [TJ For equal gap energies in the box and 
the reservoir (A r = At = A) the activation energy SE is 
determined by the effective charging energy of the CPB. 
Note that tunneling of an unpaired electron into the is- 
land shifts the net charge of the island by le. Thus, one 
can find SE of Eq. ([2]) as the energy difference at two 
values of the induced charge, N g = 1 and N g — 0, on the 
even-electron branch of the spectrum (see Fig. [1} : 



SE = E even (Ng = l)-E even (N g =0). 



(3) 



Here we assumed that subgap conductance due to the 
presence of an unpaired electron is negligible 18 . 

In order to find activation energy SE given by Eq. ([3]), 
we calculate the partition function Z(N g ) for the sys- 
tem, island and reservoir, with even number of electrons. 
For the present discussion it is convenient to calculate 
the partition function using the path integral description 
developed by Ambegaokar, Eckern and SchonJ^. In this 
formalism the quadratic in Q interaction in Eq. ([T]) is 
decoupled with the help of Hubbard-Stratonovich trans- 
formaion by introducing an auxiliary field ip (conjugate 
to the excess number of Cooper pairs on the island). 
Then, the fermion degrees of freedom are traced out, and 
around the BCS saddle-point the partition function be- 
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FIG. 1: (color online). Energy of the Cooper-pair box as a 
function of dimensionless gate voltage N g in units of e. Solid 
(red) line corresponds to even-charge state of the box, dashed 
(blue) line corresponds to the odd-charge state of the box. 
Here SE is the ground state energy difference between the 
even-charge state (no quasiparticles in the CPB), and odd- 
charge state (an unpaired electron in the CPB) at N g = 1. 
(We assume here equal gap energies in the box and reservoir, 
A r = A 6 = A.) 



comes 



Z ( N 9 ) = E e lwNam Jd^ J 



D<p(r)e 



(4) 



Here the summation over winding numbers accounts for 
the discreteness of the charge^, and the action S reads 
(h=l) 



S = 



dr 



2e 



— Ej cos <p(r) 



(5) 



dr fdr'a{r-r') ( 1-cos M')"^) 
u Jo V 



with P being the inverse temperature, /3 = 1/T. Here 
Cgeom is the geometric capacitance of the CPB which de- 
termines the bare charging energy E c = e 2 /2C gcom ; and 
Ej is Josephson coupling given by Ambegaokar-Baratoff 
relation. The last term in Eq. ([5]) accounts for single elec- 
tron tunneling with kernel a(r) decaying exponentially 
at r > A -1 [lj|- For sufficiently large capacitance the 
evolution of the phase is slow in comparison with A -1 , 
and we can simplify the last term in Eq. {5| 



(i Jo V 

3tt 2 



1 



128 2ne 2 R N A J n \ dr 



(6) 



It follows from here that virtual tunneling of electrons 



between the island and reservoir leads to the renormal- 
ization of the capacitance^ 



a 



C ■ Cgcom H~ 



3tt 1 
32 R N A' 



(7) 



Within the approximation (J6j> , the effective action ac- 
quires simple form 



Sett = dr 
h 



C ( <p(r) 
2 I 2e 



Ej cos tp{r) 



(8) 



To calculate Z(N g ) one can use the analogy between the 
present problem and that of a quantum particle moving 
in a periodic potential, and write the functional integral 
as a quantum mechanical propagator from ip^ = (pq to 
iff = ipo + 2-7rra during the (imaginary) "time" (3 



Dip(r) exp(-Seff) = (tpo\e 



(9) 



The time-independent "Shrodinger equation" corre- 
sponding to such problem has the formal 



q2 

ir, ^>(s) = E^(p), H eS =(-4:E c -^-E J cosip 



Here E c denotes renormalized charging energy 



E r . = 



1 + &g T E < 



(10) 



(11) 



32 Ht A 



One can notice that Eq. (jTUJ) corresponds to well-known 
Mathieu equation, for which eigenfunctions 'fyk.sity) are 
known™. Here quantum number s labels Bloch band 
(s = 0,1,2,...), and k corresponds to the "quasi- 
momentum" . By rewriting the propagator §§§ in terms 
of the eigenfunctions of the Shrodinger equation (JTUJ) we 
obtain 



{ip \ e^ 6 -* |<po + 27rm) 

= J2^o\k){k\, 



(12) 



\k'}{k , \(po + 2Trm) 



k,k' 



^*fe !S (^o)*fc, s (^o + 27rTO)exp(-/3£; s (fc)). 



Here E s {k) are eigenvalues of Eq. fTI 

According to the Bloch theorem, the eigenfunc- 
tions should have the form 4'fe )S (</j) = e lk,p ^ 2 Uk tS ((p) 
with Ufc, s (^) being 27r-periodic functions, Uk >s (tp) = 



Uk,s(ip + 27r). We can now rewrite Eq. as 



5E/E C 



Z(N 9 ) = E 



m— — oo 



d(po x 



x ^fe )S (^o)*fc, s (^o+27rm) exp (-(3E s (k)) 

k,s 
oo 

= E exp (iV,)). (13) 

s=0,l 

The eigenvalues E s (N g ) are given by the Mathieu char- 
acteristic functions Ma(t, q) and Mb (r, At N g — 
and N g = 1, the exact solution for the lowest band reads 

(14) 



£ o (iV s = 0) = E c M A (0-^f- J , 



£ (^ = 1) = S c Ma(1,-^ 



The activation energy SE can be calculated from Eq. (fl~3|) 
by evaluating free energy at T = 0: 



SE = E c 



M A 1 



2E r 



-M A [ 0, 



2E r 



(15) 



The plot of SE as a function of Ej/2E C is shown in 
Fig. ([2]) . Even-odd energy difference SE has the following 
asymptotes: 



E c - ^Ej, 



EJ2E C 4Z\, 



SE- r 



These asymptotes can be also obtained using perturba- 
tion theory and WKB approximation, respectively. 

As one can see from Eq. (jT5"j) . SE can be reduced by 
quantum charge fluctuations. For realistic experimen- 
tal parameters^ A « 2.5K, E c w 2K and g T w 2, we 
find that even-odd energy difference is 15% smaller 
with respect to its bare value, i.e. <5i? rs 1.45K and 
jj^bare ^ ^ 7j£ Since the reduction of the activation en- 
ergy by quantum fluctuations is much larger than the 
temperature, this effect can be observed experimentally. 
The renormalization of SE can be studied systematically 
by decreasing the gap energy A, which can be achieved 
by applying magnetic field B Q . The dependance of the 
activation energy SE on A(£?) in Eq. ([15]) enters through 
the Josephson energy E 7 , which is given by Ambegaokar- 
Baratoff relation, and renormalized charging energy E c 
of Eq. (jn]|. 

The renormalization of the discrete spectrum of charge 
states in the CPB becomes more pronounced in the 
strong tunneling regime. However, the adiabatic approx- 
imation leading to the effective action S ff (JSJ) is valid 
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FIG. 2: Dependence of the even-odd energy difference SE on 
the dimensionless parameter E j/2E c . 



when the evolution of the phase is slow, i.e. the adia- 
batic parameter uj f /A is small. (Here lo , is the plasma 
frequency of the Josephson junction, u>j ~ -J E c Ej.) 
Thus, at large conductances gx the adiabatic approxima- 
tion holds only when the geometric capacitance is large 
Cgeom ^ e 2 <? T /A. Under such conditions the renormal- 
ization effects lead to a small correction of the capaci- 
tance, see Eq. ([7]). If w ; /A > 1, the dynamics of the 
phase is described by the integral equation ([5]), and re- 
tardation effects have to be included. 

In the similar circuit corresponding to the Cooper-pair 
box qubiti^ it is possible to achieve strong tunneling 
regime g T 3> AC gcom /e 2 , and satisfy the requirements 
for adiabatic approximation (u> , / A «C 1). In this circuit 
a single Josephson junction is replaced by two junctions 
in a loop configuration^. This allows to control the ef- 
fective Josephson energy using an external flux <fr x . (For 
the CPB qubit the Josephson energy Ej in Eq. ([8]) should 
be replaced with Ej{$ x ) = 2E° cos (n$ x /$ ); here $ 
is the magnetic flux quantum, $o = h/2e, and E° is 
the Josephson coupling per junction.) In this setup even 
at large conductance g T ^> AC gcom /e 2 one can decrease 
luj ~ y/ E c Ej (J&x) by adjusting the external magnetic 
flux to satisfy wj/A <C 1. Under such conditions the 
quantum contribution to the capacitance C (see Eq. 0) 
becomes larger than the geometric one, while the dy- 
namics of the phase is described by the simple action of 
Eq. ||5J). It would be interesting to study experimentally 
the renormalization of the discrete energy spectrum of 
the qubit in this regime. We propose to measure, for ex- 
ample, the even-odd energy difference SE. In this case 
SE is determined by the conductance of the junctions 
<? T , superconducting gap A, and magnetic flux <& x , and 
is given by Eq. (15]) with E c « 32A/3g T , see Eq. (fTT|) . 
and Ej = 2E° cos (ir<I> x /<S> Q ). 

In conclusion, we studied the renormalization of the 
discrete spectrum of charge states of the Cooper-pair box 
by virtual tunneling of electrons across the junction. In 
particular, we calculated the reduction of even-odd en- 
ergy difference SE by quantum charge fluctuations. We 
showed that under certain conditions the contribution of 
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quantum charge fluctuations to the capacitance of the NSF grants DMR 02-37296, and DMR 04-39026. 
Cooper-pair box may become larger than the geometric 
one. We propose to study this effect experimentally using 
the Cooper-pair box qubit. 
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